HU-EP-11/37, ITEP-LAT/2011-07 



O 

(N 
o 

Q 



in 

(N 
> 
in 
m 
r-« 

od 
o 



Landau gauge gluon and ghost propagators at finite temperature from quenched 

lattice QCD 

R. Aouane, 1 V. Bornyakov, 2 E.-M. Ilgenfritz, 3 V. Mitrjushkin, 4 M. Miiller-Preussker, 1 and A. Sternbeck 5 

1 Humboldt- Universitat zu Berlin, Institut fur Physik, 12489 Berlin, Germany 

^Institute for High Energy Physics, 142281, Protvmo, Russia 
and Institute of Theoretical and Experimental Physics, 117259 Moscow, Russia 
3 Humboldt- Universitat zu Berlin, Institut fur Physik, 12489 Berlin, Germany 
and Joint Institute for Nuclear Research, VBLHEP, 141980 Dubna, Russia 
* Joint Institute for Nuclear Research, BLTF, 141980 Dubna, Russia 
and Institute of Theoretical and Experimental Physics, 117259 Moscow, Russia 
5 Universitat Regensburg, Institut fur Theoretische Physik, 93040 Regensburg, Germany 

(Dated: December 31, 2011) 

The behavior of the Landau gauge gluon and ghost propagators is studied in pure 5(7(3) gauge the- 
ory at non-zero temperature on the lattice. We concentrate on the momentum range [0.6, 2.0] GeV. 
For the longitudinal as well as for the transverse component of the gluon propagator we extract 
the continuum limit. We demonstrate the smallness of finite-size and Gribov-copy effects at tem- 
peratures close to the deconfinement phase transition at T = T c and within the restricted range 
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to the phase transition we propose some combinations of Dl{o) signalling the transition much like 
"order parameters". 
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I. INTRODUCTION 

It is commonly believed that hadronic matter 
at high temperature undergoes a phase transition 
into another phase, traditionally called "quark-gluon 
plasma" . At present, strong efforts are made at RHIC, 
BNL and at the LHC, CERN to establish undeniable 
experimental signatures in the final states of heavy- 
ion collisions indicating that matter had undergone 
evolution close to or beyond this transition. The ex- 
istence of such a transition has been concluded long 
time ago from Hagedorn's thermodynamical model [l| 
and later has been one of the first crucial forecasts of 
lattice QCD (LQCD). The latter uses a formulation 
of non-Abelian gauge theory which is amenable to ab- 
initio numerical non-perturbative computations. This 
formulation also opens the way for analytical calcula- 
tions at strong and weak coupling. LQCD calculations 
can provide estimates for the transition temperature, 
the equation of state close to T c and above and other 
features and experimental observables. For a recent 
review see 0]. 

In recent years, another powerful non-perturbative 
approach has been developed based on Dyson- 
Schwinger equations (DSE) [3H6| and functional 
renormalization group equations (FRGE) 0, H[ . The 
main focus was first to find a field theoretical, model- 
independent description of quark and gluon confine- 
ment in terms of the infrared behavior of gauge- 
variant Green's functions, in particular of the Landau 
or Coulomb gauge gluon and ghost propagators. This 
should allow to confirm or disprove confinement sce- 



narios as proposed by Gribov and Zwanziger [9Hll| 
and Kugo and Ojima [H, E3- Landau gauge gluon 
and ghost propagators have been intensively studied 
for zero temperature with DSE and FRGE (see, e.g., 
(l4j and citations therein). On the lattice these prop- 
agators have been computed by several groups (see 
EB] for our own recent computations and refer- 
ences to earlier work by other groups) . 

If these propagators encode confinement, they 
should also be considerend in LQCD studies at non- 
vanishing temperature (see [17H20l | and for more re- 
cent work [2l| - |24| ). Complementary to this the tem- 
perature dependence has also been studied in the 
framework of DSE [H-[28|. 

In the recent past we also have extended our lattice 
computations of the gluon and ghost propagators at 
zero temperature to the case of non-zero temperature. 
First results have been obtained without 29| and also 
with Nf — 2 dynamical fermion flavors 30(. In this 
paper we will focus on results for pure SU(3) gauge 
theory known to have a first order finite temperature 
phase transition. We concentrate on the continuum 
limit within a restricted range of momenta. For this 
range finite-size or Gribov-copy effects turn out to be 
small. Moreover, a noticeable sensitivity of the longi- 
tudinal component of the gluon propagator with re- 
spect to the deconfining phase transition is observed. 
We show that certain ratios of this component may 
serve as useful indicators ( "order parameters" ) for this 
transition, thus complementing the information ob- 
tained from the ever popular Polyakov loop. 

The paper is organized as follows. In Section ILT1 wc 
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describe the setup of our lattice Monte Carlo simula- 
tion. In Section IIIII we review the basic definitions of 
the gauge- variant propagators on the lattice, modified 
to finite temperature. In Section[TV]we present results 
for the gluon and ghost propagators for various tem- 
peratures. The signal of the phase transition is not as 
strong as one might have expected. Nevertheless, as 
it is said above, the longitudinal gluon propagator al- 
lows us to define ratios which give a clear signal at the 
deconfinement phase transition. In Sections IVl and PvTI 
we analyze the finite- volume and Gribov-copy effects, 
respectively. In Section fVIII we investigate then scal- 
ing properties for varying lattice spacing a, keeping 
the temperature and the volume fixed. This allows us 
to extrapolate our data to the continuum limit. Fi- 
nally, in Section [Villi we shall draw our conclusions. 



II. SETUP OF THE LATTICE SIMULATIONS 

We have generated SU (3) pure gauge field configu- 
rations on a four-dimensional lattice of size TV 3 x N T 
with periodic boundary conditions employing stan- 
dard Monte Carlo simulations using the Euclidean 
path integral weight ~ exp (— Sw), where Sw denotes 
the Wilson one-plaquette action 



Sw = P 
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go is the bare coupling constant and U xfi <E SU(3) de- 
notes the link variables. The imaginary-time extent 
corresponds to the inverse temperature T^ 1 — N T a, 
where a{0) is the lattice spacing. For generating the 
gauge field ensemble we have used the standard hy- 
brid over-relaxation algorithm, with a step of 4 mi- 
crocanonical over-relaxation sweeps followed by one 
heatbath step [3l|, [HJ . In both steps a decomposition 
of SU(3) link variables into SU(2) matrices, as pro- 
posed in [33]], was applied. 0(2000) combined ther- 
malization sweeps were allowed between the individ- 
ual measurements of the propagators. 

In order to determine the temperature dependence 
of the gluon and ghost propagators, in a first step 
we kept the lattice spacing fixed (and, as we shall 
see, sufficiently small) while varying N T . As a ref- 
erence value we have chosen (3 = 6.337 providing 
a ~ 0.055 fm (in accordance with [Hj]). This /3-value 
corresponds, for N T — 12, to a temperature very close 
to the temperature T c characteristic for the deconfine- 
ment phase transition in a lattice with a linear spatial 
extent N a a(f3 = 6.337) = 48 a ~ 2.64 fm. According 
to Ref. [HI it has been fixed by interpolating with the 
help of the fit formula 
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TABLE I: Temperature values, lattice size parameters, val- 
ues of the inverse bare coupling /3, the lattice spacing a in 
units of GeV and fm, the number n con f of independent 
lattice field configurations and the number n copy of gauge 
copies used throughout this study. 



where /3 C (N T , oo) corresponds to the thermodynamic 
limit and h denotes a fitted coefficient (h < 0.1). 
N a = 48 guarantees a reasonable aspect ratio over 
the whole temperature range T/T c = 12/N T e 
[12/18, 12/4] and permits to reach three-momenta be- 
low 1 GeV. 

As a second step, we decided to study systematic ef- 
fects as there are finite- volume effects (cf. Section fV|) . 
Gribov copy effects (cf. Section IVI[) and the scaling 
properties (cf. Section IVIIj) in order to extrapolate to 
the continuum limit a — ¥ for a couple of momen- 
tum values. For the two latter studies we varied a(/3) 
while having kept constant the physical spatial volume 
(2.7 fm) 3 as well as two representative temperature 
values (T ~ 0.86 T c and T ~ 1.20 T c , respectively). 

A compilation of the lattice sizes (N T x AT 3 .) and (3- 
values together with the number of independent lat- 
tice configurations generated for this study can be 
found in Table U 



III. GLUON AND GHOST PROPAGATORS 

For determining the gluon and ghost propagators 
we have to fix the gauge. Under local gauge transfor- 
mations {g x } the link variables transform as 



/3 c (iV T ,A CT ) = /3 c (iV T ,oo) 



N„ 



U Xfl & = glU xll g x+fl , 



g x E SU(3) . (1) 
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In order to satisfy the Landau gauge transversality 
condition 

Vp^ = (2) 
with the lattice gauge potentials 

A^(X + jl/2) = ^—(U XIJl - U^) \traceless (3) 

it is sufficient to maximize the gauge functional 

F u[9] = lj2^ eTl 9 x U XIM gl +jl (4) 

XyfJ, 

with respect to g x . What concerns the Gribov non- 
uniqueness problem for solutions of the gauge condi- 
tion Eq. ([2} we adopt the strategy of finding gauge 
copies being as close as possible to the global maxi- 
mum of Fjj [g] [H, H3| as already practized in [l5|, [l6l 
l38l - l4l| . This prescription has been shown to provide 
correct results for Landau gauge photon and fermion 
propag ators within compact U(l) lattice gauge theory 
|42h45| . Very efficient for this aim is the simulated 
annealing (SA) algorithm combined with subsequent 
overrelaxation (OR) iterations [H, [H, [H-Hl]. The 
SA algorithm generates gauge transformations {g x } 
randomly with a statistical weight ~ exp(Fu[g]/T sa ). 
The "temperature" T sa is a technical parameter which 
is monotonously lowered in the course of 3500 SA sim- 
ulation sweeps (actually, these are heatbath updates). 
Also, for better performance, a few microcanonical 
steps are applied after each heatbath step. In fact, 
we start with T sa — 0.45 and decrease this parameter 
down to T sa = 0.01 in equal steps after each combined 
sweep. Finally, in order to satisfy the gauge condition 
Eq. © with a local accuracy of 

max<Ke Tr[V nA X nV v At v } < e , e = 10~ 13 (5) 

X 

we employ the standard OR procedure. Except for the 
study of the influence of Gribov copies (cf. Section lVH) . 
we carry out only one such attempt per configuration 
to fix the gauge. As in our previous studies, we call 
the corresponding (first trial) gauge copy "first copy" 
(/c). 

The gluon propagator is defined in momentum 
space as 

Df v (q) = (K(k)Al(-k)), (6) 

where (■ • • ) represents the average over configurations, 
and A^(k) denotes the Fourier transform of the gauge- 
fixed gluon field ([3]) depending on the integer-valued 
lattice momentum k^ (p = 1, . . . , 4). The latter is re- 
lated to the physical momentum (for the Wilson pla- 
quette action) as 

q,(k») = - sin (^) , k, e {-N,/2, NJ2] , (7) 



where (Ni,i — 1,2,3; N4) = (N a ;N T ) characterizes 
the lattice size. 

For non-zero temperature it is convenient to split 
the propagator into two components, the transverse 
Dt ( "chromomagnetic" ) (transverse to the heatbath 
rest frame) and the longitudinal Dl one ("chromo- 
electric"), respectively, 

DfM = 5 ab (P^D T (ql f) + P^D L {ql, f )), (8) 

where qn plays the role of the Matsubara frequency, 
which will be put to zero lateron. For the Landau 
gauge, the tensor structures P^ L represent projectors 
transverse and longitudinal relative to the (p = 4)- 
direction 

p£, = (i - - M (<V - ^) - (9) 

For the propagator functions Dt,l we find 

DT= 2W g /E4 b (*)^ 8 (-*0 - ^Ai(k)AZ(-k)^ 

(11) 

and 

DL = W g { 1 + f) (^( fe K(-*)) . (i 2 ) 

where the number of generators N g = N 2 olm . — 1 for 
iVcolor = 3. The zero-momentum propagator values 
can be defined as 

^(°) = ^E^^c))- (1 3 ) 

& i—i 

Dl(0) = ±(Ai(p)AZ(p)). (14) 

Notice that - at least for large enough j3 - the Lan- 
dau gauge gluon propagator is expected to depend on 
the Z(3)-sectors into which the Polyakov loop spa- 
tial averages can fall [igj]. Therefore, before carrying 
out the SA gauge fixing procedure we always apply a 
Z(3)-flip as described in Section IVT1 but with respect 
to the 4-th direction. It ensures the phases of the 
corresponding Polyakov loop averages to fall into the 
interval (— tt/3, 7r/3]. 

The Landau gauge ghost propagator G(q) and its 
dressing function J(q) are defined as follows. 

G ab {q) = a 2 ^(e- 2 ™( fe / w H*-y)[ M -i]afc) ; 
x,y 

= 8 ab G(q) EE 5 ab J(q)/q 2 , (15) 

where q 2 ^ and (k/N)-(x-y) = £ M k^x-y^/N^. 
M denotes the lattice Faddeev-Popov operator cor- 
responding to the gauge field definition © and the 
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related gauge functional i.e., 



with 



A% = mcTr[{T a ,T b }(U x ^ + U x ^% 



,v 

;b 
'x,y 

J x,y 



B a }„ = 2 • D\eTi[T b T a U x ^}, 



C±, = 2 ■ meTr[T a T b U x . 



where T a (a = 1, . . . ,N g ) are the Hermitian genera- 
tors of the su(3) Lie algebra satisfying Tr[T Q T b ] = 
5 ab /2. In order to invert M we use the conju- 
gate gradient (CG) algorithm with plane- wave sources 
ip c with color and position components ipc( x ) = 
<5" exp (2iri(k/N) ■ x). Actually, we apply a pre- 
conditioned CG algorithm to solve the equations 
M£y4 ' (y) = V'c( x )> where as pre-conditioning ma- 
trix we use the inverse Laplacian A -1 with a color- 
diagonal substructure [Ho|, HU . 

In order to study hypercubic lattice artifacts, we 
have analyzed the influence of the choice of momenta 
on the behavior of the gluon propagator. When com- 
paring on-axis with diagonal momenta (for /3 = 6.337 
and in the lower momentum range) we found only 
small but nonetheless systematic deviations due to 
the hypercubic lattice geometry. To maximally reduce 
them the so-called cylinder cut [52] 



E fc M-i(EM 2 < 



(17) 



with A: 4 = and c — 3 has been applied to all our 
data. 



IV. RESULTS: GLUON AND GHOST 
PROPAGATORS VERSUS TEMPERATURE 

In Fig. [T]we display the multiplicatively renormal- 
ized propagators D^{q) and Dx(q) as functions of the 
three-momentum (q = \q\, q^ = 0) for j3 = 6.337, ob- 
tained with N a = 48 and different N T , i.e. for temper- 
ature values varying from T = 0.65 T c up to T ~ 3 T c . 
For details we refer to the upper section of Tablefl] The 
renormalization condition is chosen such that Dl,t 
take their tree level values at the subtraction point 
q = fi. We choose ft = 5 GcV in order to be close to 
the perturbative range and still reasonably away from 
our lattice cutoff (q m ax — 2v3/a ~ 12.4 GeV). 

One can see from Fig. Q] that the temperature de- 
pendence of both Dl and Dt becomes weaker with in- 
creasing momentum. This weakening proceeds faster 
for Dt than for Dl- The ultraviolet regions of Dt 
and Dl turn out to be "phase-insensitive". This ob- 
servation was also reported in 30]. More precisely, 
while the temperature changes from its minimal value 



to our maximal one, the change of Dt is less than 5% 
for q > 2.2 GeV, while for Dl this is guaranteed for 
q > 2.7 GeV. For T ^ T c Dl shows a comparatively 
weak temperature dependence also at small momenta. 
This changes drastically as soon as T ^ T c . In con- 
trast to that Dx{q) changes monotonously with T in 
the infrared region. This can be seen in more detail 
from Fig. [2l There we show the temperature depen- 
dence of DL{q) (left panel) as well as of Dt(<?) (right 
panel) for six selected momenta in the range up to 
1.6 GeV. 

One can see that Dl at fixed momentum shows 
strong variations in the neighbourhood of T c . It is 
rising with T below T c and sharply drops around T c . 
This behavior looks most pronounced for zero mo- 
mentum and gets progressively weaker at higher mo- 
menta. For the lowest momenta we observe maxima 
at T = 0.86 T c . It remains open, whether the maxima 
arc shifted away from the transition temperature with 
increasing volume. 1 

In any case, our data confirms that the infrared part 
of Dl (p) is strongly sensitive to the temperature phase 
transition [2^, |3(|. It may serve to construct some 
kind of order parameters characterizing the onset of 
the phase transition, as we will propose below. In 
contrast to that, Dt is ever decreasing and varying 
smoothly across T Cl showing no visible response to 
the phase transition at all momenta. 

We fit the momentum dependence in the range [0.6 : 
8.0J GeV with a Gribov-Stingl interpolation formula 
[91.I531] used in Refs. and derived lateron in the 

so-called "Refined Gribov-Zwanziger" approach (55l . 

M 



D(q) = 



c (1 + dq 2n ) 
(q 2 + r 2 ) 2 + b 2 



(18) 



Expected logarithmic corrections needed for the ultra- 
violet limit have been neglected here (for a thorough 
discussion see (H)])- We put throughout n = 1. In 
a first attempt we have left b varying. We obtained 
values compatible with 6 = except for D^iq) at the 
highest three temperature values inspected. There- 
fore, in all other cases we have repeated the fits with 
fixed 6 = and obtained x^-values reasonably below 
2.0. The fit parameters can be found in Table El 2 
Since we expect to see a plateau and even a bend 
over for Dt at momenta below our minimal ones the 
parameter b might become nonzero also at lower tem- 
peratures. This would then correspond to a complex 
effective mass parameter. 



1 For SU(2) gauge theory the maximum of D^(0) was recently 
reported 12411 to move away from the transition with decreas- 
ing lattice spacing. 

2 Note that for 6 = Eq. d 1 S I) is equivalent to the interpolation 
formula D(g) = {qi ~l si) + {q i+ 5 iyi ■ 
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FIG. 1: Temperature dependence of the longitudinal (l.h.s.) and the transverse (r.h.s.) gluon propagator for /3 — 6.337 
and a spatial lattice size N a = 48. 
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(l.h.s.) and the transverse one, Dt, (r.h.s.) 
k 3 ,ki). /3 = 6.337 and N a = 48. 



vs. temperature for a few low 



We have tried to form quantities constructed from 
the gluon propagator which can serve as indicators for 
the deconfincmcnt transition. First, we plot the ratio 

X =[D L (0,T)-D L (q,T)}/D L (0,T) (19) 

as a function of T/T c in the left panel of Fig. [3J We 
observe that all the curves labelled by the momen- 
tum 4-tuples in the legend show approximate plateaux 
below T c . Then, passing the phase transition they 
suddenly fall off with slopes becoming slightly smaller 
with increasing momentum, but still with visible tem- 
perature sensitivity. This means that x can be used 
as an indicator for the deconfinement transition and, 
moreover, the transition can be traced even at rather 
high momentum. This was not so clear from the 
l.h.s. of Fig. [2 where the behavior of Dl at higher 
momenta looks rather smooth. 

From the behavior of x^ & t least in the interval 
0.65 T C ^ST ^ST c and at low momentum, one can con- 
jecture the factorization 

D L (q;T)~A(q)-B(T). (20) 



Then, as long as the temperature T varies in the given 
interval, the change of Dl can be described by a mo- 
mentum independent rescaling. This is a rather non- 
trivial property from which further conclusions can 
be drawn. For example, in the interpolation formula 
(fT8|) above, we should find the mass parameter r 2 and 
the parameter d to be (approximately) temperature 
independent as long as T < T c . 

From the left panel of Table [XT] one can see that 
this is true for the parameter d which varies within 
error bars. The variation of parameter r 2 is up to 
30%. This comparatively large variation might be ex- 
plained by the fact that the propagators were fitted 
over a wide range of momenta, specifically from 0.6 to 
8 GeV, while the factorization we expect to hold only 
at low momenta. Indeed, a fit in the range up to 2.5 
GeV which includes also the zero momentum value 
has shown r 2 to become approximately constant. It 
remains to be seen, whether this behavior in the in- 
frared region survives the continuum and the thermo- 
dynamic limites which goes beyond the scope of this 
paper. 
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Parameters 
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1.398(62) 


0.133(8) 


6.15(34) 


0.93 


2.97 


4 


5.327(159) 


0.0 


0.045(2) 


17.15(73) 


0.51 




0.927(126) 


2.559(33) 


0.100(6) 


7.58(41) 


1.01 



TABLE II: Results from fits with Eq. (fT8)) (n = 1) for Dl (l.h.s.) and Dt (r.h.s.) corresponding to the Monte Carlo 
data shown in Fig. [T] (/3 = 6.337, N a = 48). The fit range is [0.6 : 8.0] GeV. The values in parentheses provide the fit 
errors. The boldface printed 6-values indicate that they are fixed to zero. 
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Let us consider another ratio 



D L (0,T)-D L (q,T) 
£>l(0, T mi „) — D^{q,T min ) 



Train — 0.65 T c , 

(21) 

which according to the factorization (|2"0|) should be 
approximately momentum independent. Indeed, this 
can be seen from the right panel of Fig. [3] Moreover, 
a(q, T) should resemble qualitatively the temperature 
dependence of Dl at q = 0. Close to T c , however, 
a falls off reaching very small values at higher tem- 
peratures (around 2T C ). Therefore, we conclude that 
both quantities \ (ceasing to be constant) and a (with 
its strong fall off) signal the finite-temperature transi- 
tion. It remains to be seen, whether they also map out 
the (pseudo)critical behavior in unquenched QCD. 

Let us note that our volumes are not large enough to 
study the infrared asymptotic behavior. Moreover, at 
the lowest momenta we expect systematic deviations 
due to finite-size effects, lattice artifacts, and Gribov 
copy effects. This concerns also the parameters \ an d 
a because of their dependence on the value D^(q = 
0). The systematic effects will be discussed to some 
extent in Sections [V] to I VIII in order to identify the 



momentum range, where they play only a negligible 
role. 

Summarizing, in agreement with findings in other 
recent investigations [2lT - |24l . l30l ] we observe the 
strongest response to the phase transition to occur 
in the gluonic chromoelectric sector (the longitudinal 
propagator) rather than in the gluonic chromomag- 
netic one (the transverse propagator). 

We have also computed the ghost propagator ac- 
cording to Eq. (|15p , restricting it for simplicity to the 
diagonal three-momenta and vanishing Matsubara fre- 
quency, k u = (fc, k, k, 0) with k = 1, . . . , 7. The data 
are again normalized at fj, — 5 GeV, such that the 
ghost dressing function equals unity at q = fi. The 
result for the latter function is displayed in Fig. H) 
In comparison with the gluon propagator we see the 
ghost propagator to change relatively weakly with the 
temperature. 3 This is in agreement with the observa- 
tion in [21] . An increase becomes visible at tempera- 



3 Note the use of a linear scale at the vertical axis in Fig. [4^, 
in contrast to the logarithmic one in Fig. [T] 



7 



ture values T > 1AT C for the lowest momenta studied 
(see Fig.|4jD). The relative insensitivity with respect to 
the temperature is the reason why we will not further 
consider the ghost propagator in what follows. 



V. FINITE- VOLUME EFFECTS 

In order to estimate finite-volume effects we com- 
pare the data shown before with data obtained on 
even larger spatial volumes while keeping fixed the 
coupling (at j3 = 6.337) and two temperature values, 
T = 0.86 T c (confinement) and T = 1.2T C (deconfine- 
ment), respectively. The linear spatial extent varies 
from 48a = 2.64 fm to 64a = 3.52 fm (see also the 
middle section in Table |T| . 

In Figs. [5]and[H]we show the corresponding plots for 
Dl and Dt, respectively. In all four cases we observe 
the effects to be small for momenta above 0.6 GeV. 4 
For lower momenta, especially at zero momentum, 
systematic deviations become more visible. With in- 
creasing volume the infrared values of Dl seem to rise, 
whereas for Dt the opposite is the case. This behavior 
has already been reported for pure gauge theories in 
[H[52l for SU(2) and in [H for SU{3), respectively. 



VI. GRIBOV COPY EFFECTS 

In order to study Gribov copy effects we compare 
"first" , i.e. randomly occuring copies (fc) with "best" 
copies (6c). The latter were produced as follows. 

We searched for copies within all 3 3 = 27 Z(3) sec- 
tors characterized by the phase of the spatial Polyakov 
loops, i.e. Polyakov loops in one of the three spatial 
directions. For this purpose the Z(3) flipping oper- 
ations [10, Eo[ were carried out on all link variables 
U Xt i (i = 1,2,3) attached and orthogonal to a 3d 
hyperplane with fixed Xi by multiplying them with 
exp {±2-7ri/3}. Such global flips are equivalent to non- 
periodic gauge transformations and do not change the 
pure gauge action. For the 4th direction, we stick to 
the sector with \arg P\ < 7r/3 which provides maximal 
values of the functional dl]) at the /3- values considered 
in this section [3(J. Thus, the flip operations com- 
bine for each lattice field configuration the 27 distinct 
gauge orbits of strictly periodic gauge transformations 
into one larger gauge orbit. 

The number of copies actually considered in each 
of the 27 sectors depends on the rate of convergence 
(with increasing number of investigated copies) of the 



4 Below T c the transverse propagator changes by less than 
12 %, the longitudinal one by less than 5 %. Above T c the 
transverse propagator varies by less than 8 % and the longi- 
tudinal one by less than 11 %. 



propagator values assigned to the best copy, in par- 
ticular at zero momentum. From our experience with 
SU (3) theory [30| we expect that the effect of consid- 
ering gauge copies in different flip-sectors is more im- 
portant than probing additional gauge copies in each 
sector. For this reason and to save CPU time we 
have considered one gauge copy for every Z(3)-sector; 
therefore, in total n copy = 27 gauge copies for every 
configuration. 

To each copy the simulated annealing algorithm 
with consecutive overrelaxation was applied in order 
to fix the gauge. We take the copy with maximal 
value of the functional (UJ as our best realization of 
the global maximum and denote it as best ("6c") copy. 

The parameters of the SA algorithm in the study 
of Gribov copies were slightly different from those de- 
scribed above in Section ITTT1 2000 SA combined sim- 
ulation sweeps with a ratio 11:1 between microcanon- 
ical and heat bath sweeps were applied starting with 
T sa = 0.5 and ending at T sa = 0.0033. 

Since this procedure is quite CPU time consuming 
we restricted this investigation to coarser lattices 6 x 
28 3 and 8 x 28 3 with larger lattice spacing, such that 
both the temperature values T — 0.86 T c and T = 
1.20 T c , respectively, as well as the physical 3d volume 
(2.64 fm) 3 were approximately reproduced. 

In Figs. [7] and [8] we compare 6c with fc results for 
the gluon propagators Dl and Dt, respectively. As 
one can see, Dl is almost insensitive to the choice of 
Gribov copies (at least, for the comparatively small 
values of Nj- we consider), as has been already re- 
ported in [29[ for the SU(2) case and in [3(| for the 
SU(3) case. On the contrary, the transverse propa- 
gator is strongly affected in the infrared region. This 
observation is independent of the temperature. More- 
over, we see that the transverse gluon propagator val- 
ues in the infrared become lowered for 6c compared 
with fc results. These observations resemble those 
made already in [2!| and [30| . 

The main conclusion of this section is that Gribov 
copy effects may be neglected for all nonzero momenta 
in the case of the longitudial propagator (at least, for 
comparatively small values of N r ), and for momenta 
above 800 MeV in the case of the transverse propaga- 
tor. The momentum range where the last statement 
is true might depend on the temperature. 



VII. SCALING AND CONTINUUM LIMIT 

In order to check for good scaling properties we have 
used the same reference values for the temperature 
below and above T c as discussed before (i.e., 0.86 T c 
and 1.20 T c ). We kept also the spatial volume fixed at 
(2.7 fm) 3 . We compare the renormalized propagators 
at four different values for the lattice spacing a(f3) (see 
Table [XJ) . Our results are displayed for the momentum 
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range up to 1.5 GeV in Fig. [9] for Dl and in Fig. [TU] 
for Dt, respectively. Gauge fixing has been carried 
out as originally described in Section IIIII 
We provide the renormalization factors for 

D L>T (q, ft) = Z L>T (a, fi) D L ar T e (q, a) (22) 

in the left panel of Table Hill As expected the Z-factors 
of Dl and Dt approximately agree. 

From Figs. [UJ and QJJ we see that the scaling vio- 
lations happen to be reasonably small for momenta 
above 0.8 GeV. This shows that our choice of a = 
0.055 fm for f3 = 6.337 was already close to the con- 
tinuum limit. 

In order to study the a-dependence at five particu- 
lar physical momenta p we need interpolations of the 
momentum dependence in between the data points. 
For the fit within the interval 0.6 GeV < q < 3.0 GeV 
we have used again Eq. (fTB"]) with parameter b fixed to 
zero. The values of the fit parameters are displayed 
in the right hand panels of Table IIIII In all cases we 
find x 2 -values per degree of freedom around or below 
unity. 

The propagators, now interpolated to the set of se- 
lected momentum values, are shown in Fig. Qj] for Dl 
and in Fig. [T^] for Dt, respectively, as functions of the 
lattice spacing a. We show them together with the 
respective fit curves 

D(a;p) = D + B ■ a 2 (23) 

assuming only 0(a 2 ) lattice artifacts. The corre- 
sponding fit results are collected in Table IIVI The re- 
spective fit parameters Dq represent the continuum 
limit values of the propagators at the preselected mo- 
menta. 

Our lattice propagator data obtained for f3 = 6.337 
as discussed in Section IIVI can now be compared with 



the values extrapolated to the continuum limit. This 
is shown in Fig. 1131 In more detail, we can compare 
the continuum extrapolated values at some lower mo- 
mentum - say at q = 0.70 GeV - with those obtained 
from a{j3 = 6.337) = 0.055 fm and interpolated to 
the same momentum. Then we find deviations be- 
ing smaller than 4%. Thus, we are really justified 
to say that the results obtained for (3 — 6.337 in the 
given momentum range are already very close to the 
continuum limit. The continuum limit extrapolated 
propagators can be easily interpolated with formula 
Eq. (|T5| . The results are drawn in Fig. [14J 

We conclude that for the higher /3-values and the 
momentum range considered in this paper we are close 
to the continuum limit. Moreover, systematic effects 
as there are finite-volume and Gribov copy effects 
seem to be negligible for momenta above 0.8 GeV. 



VIII. CONCLUSIONS 

We have presented lattice results for the Landau 
gauge gluon and ghost propagators computed in pure 
gauge SU(3) lattice theory at non-zero temperatures. 

Overall, our results agree with those published in 
[231 ] and there are hardly any deviations which are not 
due to the lattice discretization or the finite volume. 
However, our aim here was to go a step further and 
to provide results in the continuum limit for tempera- 
tures below and above the deconfinement phase transi- 
tion, and this with negligible systematic finite- volume 
and Gribov copy effects. For this to become feasible 
we had to restrict the analysis to a well-defined mo- 
mentum range around 1 GeV. 

The systematic effects were studied at two temper- 
atures, T = 0.86T C and T = 1.20T C , such that we can- 
not really tell, what happens very close to T c . Since 
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Parameters 


Z-factors 




D L fits 




D T fits 
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P 


iV CT 




Zt 


Z L 
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c(GeV 2 ) 


X% 
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c(GeV 2 ) 
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4.96(13) 
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TABLE III: Left panel: Renormalization factors Zl,t of the renormalized propagators Dt,l(?, /i) according to Eq. (|22ll . 
The renormalization point is fi = 5 GeV. Right panels: Fit parameters and Xdf f° r fits of Dl (l.h.s.) and Dt 
(r.h.s.) using the generic fit function D(q 2 ) acc. to Eq. (|18[) . but with 6 = 0. The fit range is restricted to [0.6 : 3.0] GeV. 
The fit errors are indicated in parentheses. 



for our reference value f3 = 6.337 the critical temper- 
ature T c is reached with N T — 12, there is hope that 
also in this case we keep close to the continuum limit 
and the other systematic effects are under control. 

Our results and their parametrization can be fur- 
ther used to compare with the outcome of Dyson- 
Schwingcr or functional renormalization group equa- 



tions for the gluon propagators or employed as reliable 
input in Dyson-Schwinger studies of the quark prop- 
agator. It is well-known that the non-perturbative 
continuum approaches rely on a truncated tower of 
equations for the propagators and vertex functions. 
The way of how it is truncated has a strong influence 
on the behavior, especially at intermediate momenta. 
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Parameters 


D L fits 


D T fits 


T/T c 


p{GeV) 
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D (GeV~ 2 ) 
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Do (GeV~ 2 ) 
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TABLE IV: Results of the fits for Dl (l.h.s.) and Dt (r.h.s.) as a function of the lattice spacing a using the fit function 
D(a;p) acc. to Eq. (|23[) . The errors of the fit parameters are given in parentheses. Xdf m au cases is close or well below 
unity. See also Figs. llll and [T2l 



Lattice results from first principles as those presented 
here can help to tune the truncation correspondingly. 

Concentrating on this aim we have been forced to 
choose the lattice spacing a and the linear spatial lat- 
tice extent N a such that we were prevented from going 



far towards the infrared limit. Therefore, we were not 
able to clarify the question, what the correct behavior 
is in the far infrared region. Concerning this region 
we know, that the Gribov problem is serious and still 
not completely understood. This is also the reason, 
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why we did not try in this paper to give estimates for 
the color-electric and -magnetic screening masses. 

Our results confirm that, contrary to the transverse 
gluon propagator D T and to the ghost propagator G, 
the longitudinal gluon propagator Dl is sensitive to 
the deconfinement transition. However, despite of the 
fact that we are faced with a first order phase transi- 
tion the response to it occurs relatively smooth. We 
were able to propose two parameters built from Dl 
which can be employed as indicators (kind of "order 
parameters" ) for the transition itself. One might hope 
that they will help to shed additional light on the 
transition region also in full QCD. There, at least for 
Nf — 2 quark flavours within the range of intermedi- 
ate pion masses a rather smooth crossover is expected 
(see e.g. Q and references therein, as well as [58l.[59j). 
We shall come back to this question in the near future. 
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FIG. 11: Dl vs. lattice spacing a for a set of different preselected momenta p. l.h.s. T = 0.86 T c ; r.h.s. T = 1.20 T c . 
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FIG. 12: Same as in Fig.[TT]but for D T . l.h.s. T = 0.86 T c ; r.h.s. T = 1.20 T c . 
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FIG. 13: Comparison of the renormalized propagators (l.h.s.) and Dr(q) (r.h.s.) obtained from the Monte Carlo 

simulation at /3 — 6.337 with some continuum limit extrapolated values. 
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